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Abstract 

We study the inviscid limit problem for the incompressible Navier- 
Stokes equation on a half-plane with a Navier boundary condition de- 
pending on the viscosity. On one hand, we prove the convergence of 
Leray solutions to the solution of the Euler equation. On the other hand, 
we show the nonlinear instability of WKB expansions in the stronger L°° 
and H" (s > 1) norms. 

Key words: inviscid limit problem, nonlinear instability, Navier bound- 
ary condition and boundary layers, shear flows 

1 Introduction 

1.1 The inviscid limit problem 

We consider the two-dimensional incompressible Navier-Stokes equation on the 
half-plane il = M x R+* with viscosity e > 0: 

dtu^ + • Vu'' - £/\u^ -f Vp'' ==0 in R+ X 
(NSe) : { div = in M+ X r2 

u^\t=o = in 

u"^ — {u\,U2) is the two-dimensional fluid speed and the internal pressure. 
We add two boundary conditions: first, the standard non-penetration condition 
for the component of normal to the boundary 

where n = (0, —1), and a Navier condition that depends on the viscosity e and 
that describes the tangential part of the fluid's speed on the boundary: 

{Su^ ■ n + a'u^)^^^ = on R+ x 9f7 

where G R, and Su"^ = ^(Vu*^ Vu*^) is called the strain tensor. 

Unlike the homogeneous Dirichlet (no-slip) boundary condition, u|y=o — 0; 
the Navier condition allows the fluid to slide along the boundary. Physically, 
the Navier condition is more accurate than the no-slip condition when one takes 
into account interaction at the boundary. The slip condition is therefore used 
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to homogenise the no-shp condition on rough or porous walls (see [TB] and [H]), 
and to model, for example, the flow of blood in capillary vessels that are a few 
microns wide [5S]. Mathematically, it is shown in [35] that the Navier condi- 
tion can be derived by taking the limit when the mean-free path goes to zero 
of renormalised solutions of the Boltzmann equation with Maxwell reflection 
boundary condition. 

This paper will deal with the inviscid limit problem for (NSe), ie we study 
the behaviour of a family of Leray solutions of (NSe) with Navier boundary 
condition, (w^)o<e<Eoj relative to u, a solution of the incompressible Euler equa- 
tion 

' dtv + vVv + Vq = 
(E): divz; = 

U2\y=G = 

when the viscosity e goes to 0. 

The inviscid limit problem is well understood in the whole space or with 
periodic boundary conditions [20,. However, in domains with boundaries, the 
lack of compactness due to the presence of boundary layers makes the problem 
considerably more difficult. In recent years, much progress has been made in the 
case of non-characteristic boundaries. For incompressible fluids, these appear in 
the case of injection or succion boundary conditions (see [Ml)- For related results 
for compressible models and general hyperbolic-parabolic systems, we refer to 
[23j . In this situation the boundary layer is of size e and has an amplitude C(l), 
and sharp stability and instability conditions can be shown (see for example [6] 
and for the study of the stability or instability of the Ekman layer in rotating 
fluids). 

On the half-plane, the inviscid limit problem remains mainly unsolved in 
the case of Navier-Stokes equations with the homogeneous Dirichlct boundary 
condition. Formally, one expects the boundary layer to be of size and hence 
to write a solution as 

u'' {t,x,y) ^ u\t,x,y) + u\t,x,e-^'^y) (1) 
in which solves (E) and solves the Prandtl equation 

[dtu^ + ■ Vu')i |y=o in ^ 


A major difficulty resides in the existence of such an expansion, because, while 
the Prandtl system is well-posed with monotonous initial conditions in the strip 
[0,L] X K"*" (a fact known since the 60s, [24]), it has recently been proved by 
D. Gerard- Varet and E.Dormy in [TT] that the Prandtl equation is ill-posed in 
Sobolev spaces on T x M+. This shows that the ansatz (|T]) is only formal in 
general, and E.Grenier showed in [13| that, even when such an ansatz can be 
justified, eg in analytic framework (see [32]), it may not be valid in i?^, and in- 
stability occurs on the derivatives in L°° . The instability phenomenon is linked 
to the linear instability of shear fiows for the Euler equation. 



dtui + uidxUi + U2dyUi - dyyUi = 

div u = 

U\y = = 
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For the Navier-Stokes equation with Navier boundary condition, when 
is a fixed number independent of e, the inviscid hmit problem is solved in 
framework, see [I], [1], [TT], [2]. Moreover, asympotic expansions of the form 

X, y) = u'{t, X, y) + ^/^u\t, x, e~^''^y) 

with the amplitude of the boundary layer much smaller than in ([T|), are rigor- 
ously justified in |15j and uniform conormal estimates in agreement with this 
formal behaviour are obtained in [21] . We also refer to [2] , [35] , and [3] for the 
study of special cases in 3D. 

Here, we shall be interested in the case where the Navier condition is as 
above with a(e) = ae^^/^, a ^ constant. Denoting the two components of the 
fiuid's speed by u\ and ul, the boundary conditions translate into the following: 

(NP) : w|(t,x,0) =0 
(NCae) : \dyUl{t,x,Q) = ^u\{t,x,Q) 

Our motivation for this is to understand the transition between the unstable 
Dirichlet case and the stable Navier case studied in [TS]. Note that for e small, 
our slip condition seems to approximate the no-slip Dirichlet condition, and we 
expect the formal asymptotic behaviour of the solution to replicate the ansatz 
of the Dirichlet case, ie 

u^{t, X, y) = u'{t, X, y) + u\t, x, e-^'^y) 

and hence we expect to observe some instability. 

As of now, the term "solution of the Navier-Stokes equation (NSae)" will 
designate a Leray solution of (NSe) satisfying the boundary conditions (NP) 
and (NCae) (the full definition is given in section 2). We obtain two results: 
firstly, convergence in of every sequence of Leray solutions of (NSae) to a 
smooth solution of (E) when e goes to zero, and secondly, we prove that there 
exist boundary layer profiles such that the WKB expansion, which is also of the 
form ^ in our case, is unstable in L°° . 

1.2 The L2 stability result 

We show that converges to v in L^(rj), extending the result obtained by 
D.Iftimie and G.Planas in [Tl] for a positive constant Navier boundary condition, 
ie [Su ■ n + au)tan = with a > 0. 

Theorem 1. Let be a solution of the Navier-Stokes equation (NSae), and v 
be the solution of the Euler equation (E). We assume that the initial condition 
Vq is in H'^(fl) for some s > 2. If Uq converges to vq in L'^{n) as e goes to 0, 
then for every T > 0, 

sup |l«-(t)-^;(t)||^ f^O 

te[o,T] 

In the case of a positive boundary condition, ie (NCae) with a > 0, we merely 
check that their calculations, based on energy estimates, allow us to conclude. 
The negative case (a < 0) will require little more work. 
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1.3 The nonlinear instability result 



Proving our nonlinear instability result below will closely follow an argument 
by E.Grenier, with which he proved that linear instability implies nonlinear 
instability for the Euler equation in |13) . and for some viscous boundary layers 
in [5] with B.Desjardins. We will therefore require a form of linear instability 
as a starting point: we will start with a shear flow that is linearly unstable for 
the Euler equation. 

Definition. A smooth 2D vector field Us is a shear flow if it is written as 
Us — {us{y),0) (the vector field and its first component are indifferently called 
Us). Note that a shear flow is automatically a stationary solution of the Euler 
equation. The shear flow Ug is linearly unstable if there exist /c £ M, A £ C with 



5R(A) > 0, and a function ^ G Hq 
u{t,x,y) = 



such that 



Xt^ikx 



e e 



*'(y) 



is a solution of the Euler equation linearised around Us, 

(EL) : dtu + Us-Vu + u- Vws = in M+ x 
with (NP) condition on the boundary. 

Theorem 2. Let Ug be a linearly unstable shear flow for the Euler equation. 
Setting Y — e^^^'^y, we generate a time- dependent boundary layer Tt^ as solution 
of the heat equation 



dtulit, Y) 



dlyu;{t,Y) 

u7(o,r) 
\dYu;{t, 0) 





Us{Y) 

a£~^/'^u;{t,Q) 



the boundary condition being dictated by (NCae); u^(t,y /y/e) therefore solves 
the Navier-Stokes system. Then, for any large n G N, there exist 5q and > Q 
such that for any < £ < eo we have u'^ a Leray solution of (NSae) such that 



< Ce" 



for some s > 0, and at a time T'^ ^ nln(e 



-r^ rr.e y 



u%T',x,y)-u;[T',^ 



(2) 



L°°{n) 



Moreover, for s > 1, 



u'{T',x,y) 



T^JL 



This is a nonlinear instability result in L°° in the usual sense for WKB ex- 
pansions of the form ([1}: we show that when the boundary layer profile u'^ is 
linearly unstable for the Euler equation, then the WKB expansion is unstable, 
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in the sense that arbitrarily small perturbations yield instantaneous amplifica- 
tion on [0,T^] to reach the amplitude 0(1). Note that converges to as e 
goes to 0. An explicit example of unstable boundary layer which matches the 
assumptions of Theorem [2] is given at the end of the paper. 

Our result differs from the result obtained by E.Grenier in jl3j on the in- 
viscid limit problem for Navier- Stokes equations with homogeneous Dirichlet 
boundary conditions: in one obtains SqE^^^ instead of Sq, and therefore, 
the result with Dirichlet conditions is not a standard instability result. This 
is due to the Prandtl boundary layers, the amplitude of which is the same as 
that of the internal layer. Indeed, after rescaling (t, x, y) ^ £~^l'^{t, x, y) (these 
new coordinates are indifferently denoted {t,x,y)), a solution of (NS-y/e) with 
homogeneous Dirichlet boundary conditions can formally be written as ([T]). In 
our case of Navier condition depending on ^/e, the rescaled boundary condition 
no longer depends on the viscosity, so, as shown by D.Iftimie and F.Sueur in 
[T5] . we can write 

u^{t, X, y) u;{y) + u\t, x, y) + e^/^w^ [t, x, -^j^ 

and the factor e^/^ in front of the boundary layer in the ansatz neutralises 
the one that appears when differentiating with respect to y. In the Dirichlet 
case, this compensation does not take place, and considering times that are 
O ((n — -j) ln(£^^)) leads to the apparent loss of instability. 

The proof of Theorem [2] will follow Grenier's approach, which mainly relies 
on the construction of a WKB approximate solution to (NSae) starting with 
a linearly unstable solution to the Euler equation, and the use of a resolvent 
estimate on the linearised Euler equation (our Theorem [3]). The strength of 
the method comes from the fact that the resolvent estimate only needs to be 
shown on a subspace - typically spaces of functions with one or a finite number 
or Fourier modes. This has allowed it to be used in other circumstances than 
the study of boundary layers: F.Rousset and N.Tzvetkov used it to prove trans- 
verse nonlinear instability of solitary waves for KdV and water-wave equations 
(see [31| and [30^). 

Organisation of the paper. In the next section, we prove the LP' conver- 
gence result. Theorem [1] In section 3, we first prove the resolvent estimate 
before proving Theorem [2l and finally we give an example of shear flow that 
fits the instability result. 

2 Proof of the stability (Theorem [1]) 

We remind the reader of the notion of Leray solution to the Navier-Stokes 
system: 

Definition, ■u'^ : R+ x J7 is a Leray solution of (NSae) if: 

1. u" e C^(R+,Ll) n L2([0,r],iJi) for every T > 0, where, if E{n) is a 
functional space on 17, designates the space of divergence-free vector 
fields, tangent to the boundary and belonging to E{D,). 
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2. u"^ is a weak solution to (NSae) in the following sense: we have 

- I I ■dtip + 2ay/e [ [ ■ ip + 2e [ [ Su"" : Sip 
Jm+ Jn Jm+ Jon Jm+ Jn 

- f f {u'-V^)-u'^ f u%0)-^{0) (3) 
Jr+ Jn Jn 

for every ip G 77^(R+,i?^), where A : B = X^^i.j^i.j contracted 
product of matrices A and B. 

3. For every t>Q, satisfies the following energy estimate: 

Jo Jon Jo Jn 

For a and e fixed, such solutions are known to exist and are global in time. 
Moreover, in 2D, we have uniqueness (the proof uses the Galerkin method, 
see |15)). Likewise, the existence and uniqueness of a classical global-in-time 
solution V to the Euler equation (E) when vq G H'^{n), s > 2, are well-known, 
and for every T > 0, v e L°°([0, T], iJ'*(r2)) (see [2D] and references therein). 

Integration by parts leads to the following: let f,g£ H^, with / satisfying 
(NCae). Then 

- / Af-g^^ [ f-g + 2 [ Sf:Sg (5) 
Jn Jan Jn 

This allows us to confirm that the weak formulation ([3]) contains the Navier- 
Stokes equation, the initial condition and the Navier boundary condition ((NP) 
and the divergence-free conditions being given by the choice of the spaces in 
point 1 of the definition). 

Let w"^ = if ~ V. w"^ solves the equation 

dtw' + u" ■ Vw^ -Vv- eAw" + Vip" - g) = (6) 
For i > 0, satisfies the inequality 

7;\\w'{t)\\l,,^.+ f j {w'-Vv)-w'+2a^ [ I u'-w' + 2e f f Su' : Sw' 
2 Jo Jn Jo Jan Jo Jn 

<^ll^^(0)lli.(o) (7) 

which can be seen as Q multiplied by , integrated by parts in space and 
integrated in time. Note that this is possible only because we are working in 
2D, which gives w'^ the right regularity to be used as a test function. In 3D or 
higher, ([7]) also holds; to prove it, one starts by writing the energy inequality 
then get three energy equalities: (I) by multiplying the Euler equation by v, 
(II) by again multiplying (E) by u^, and (III) by multiplying the Navier-Stokes 
equation by v. One then gets ([7]) by doing (|3])-|-(I)-(II)-(III) and using ([5]). 
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Let g > 0. For every x, y G R^, let z = .t — y. We have 



2(2: • z) = 2 



z + 



2 1 



(8) 



and the same goes for the contracted matrix product. So, from ([7]) we get 



Wilis + 4a Ve 



Jan 



-4e 



Jn 



S 



< \\w%0)\\ 



L2 ' 



Jf2 



Ve f f \v\' + e f [ \Sv\' 
Jo Jon Jo Jn 



(9) 

The left side is greater than ||iy'^(<)||22, and we estimate each term on the right- 
hand side as foUows: 

w^\\i^2 because v{t) G H'^^fl) for every t > 
and for a certain s > 2, so v{t) e C^iTi) n H^{Q), thus Vw(t) e L°°{Q); 

" /an 1"^!^ — ^7 ll'^llffi ^ certain > 0; 
• \\Sv\\'l2 < \\v\\h^ given that Sv = i(Vw +* Vu); 
so © becomes 



w^^wiii. < ii^«^(o)ii 



L2 



w||^i dT + 2 I \\Wv 



w 



LL2 



dr 



to which we apply Gronwall's lemma, and we obtain 



+ 1 / dr 



exp 2 / ||Vw 



dr 



(10) 



Fix T > 0, and, if {E{Vt), \\-\\^) is a Banach space of functions on 51, we de- 
note WfW^^E = sup4g[o,T] II/WIIb- Using the facts that v e i°°([0, T], 
and that for every t € [0,T], v{t) £ H^{D,) C Lip(f2), we rewrite as 

Ik^ WIIl^ < [lk'(0)|li. + V~e {aC^ + V^)T MI^^,] exp (2T ||«||^^Lip) 

(11) 

The right-hand side of this last inequality is bounded uniformly in t, and goes 



to as e — > 0, so we have ||w^ 



loOTi^(O) 



0. □ 



Let g < 0. Following the method above leads to a situation in which the 
coefficient in front of the exponential term after applying Gronwall's lemma no 



longer converges to zero. Instead of ([8]), we use x 



|zp -I- y • z in Q to get 



Jn 



-2e 



Sv: Sw' + 2\a\y/^ / |^o"p^-2|a|^/i 
Jo Jan 



Jan 



Regrouping all but the last two terms on the right-hand side, we write 

i Mmi. + 2s wsw'wi, <T,+T2+n 



(12) 
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/oIlv^lL. 

2 



lL2 



-eJ^i\\Siv)\\l + \\Siwn\ 



T2 < 2|a|Y^C^/Q II^'^II'hi/2(si) by trace continuity. We then use, in that 

order, the Sobolev interpolation inequahty, ||ti''^||^i/2 < ||''^'^|Il2 ll'^'^llffi j 
Young's inequahty with a positive parameter 77 > depending on e to be 
chosen later, and Korn's inequality to get 

T2 < 2\a\V~e^ f WSw^l, + {r, + 2)C7^|a|V? f \\w^l, := T' 
V Jo Jo 

We recall Korn's inequality: if / is a vector field on the half-space ft, 
tangent to the boundary, then ||V/||^2({2) < 2 ||5/||^2(f^)- If / G H^, we 
in fact have an equality. 



rt 11 ||2 
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. and finally Tg < \a\y^ J,] Jg^[\v\^ + \w^\^] < \a\y^C, /J i,.,,^! 
So becomes 

1 



Irp/ . 

2^ ' 



\\w%t)\\l.+2e 



\Sw' 



L2 



< 



\\w%0)\\i., + V~e{V^+\a\C,) / \\v\\h. 



lioo +-(r/ + 2)C^|a|Vi 



£ 1 



3\a\C\ 



\Sw' 



L2 



We fix e, and choose 77 — so that the terms with ||S'w^|j^2 cancel each 

other out. Then we take T > 0, use Gronwall's lemma and estimate the right- 
hand side independently of t e [0, T] and get 



L2 



< 



\w^0)\\l,+2V^iV^+\a\C^)T\\v\\ 



T \\v\\ 



X exp 

Again, we can conclude that ||w'^||^^^2 



0. □ 



Remark: one can study a more general Navier boundary condition depending 
on the viscosity: dyu\{t,x,0) = ae~'^ul{t, x,0). 

is then replaced by e^^^ in ([3]), (0]) and ([7]). When a > 0, in the same 
way as we got (jlip . we obtain 

\\w%t)\\^., < [\\w%0)\\^.+aAT^,s'-^]eMB^) 
When a < 0, we use ([S]), the Sobolev, Young and Korn inequalities to get 

Wn^Hml. < [\\w^O)\\l. + \a\A'^^,e'-'^]eMBrJ 

where Ar^e, A'j, ^ ^ 1 + , 9 > 0. The difference in powers of e is due to the use 
of Young's inequality when a < 0: to cancel out the terms involving ||S'w^||^2, 
one needs a coefficient rj ^ . 

Thus, in both cases, we have ranges of /3 for which the stability result 
holds: for a > 0, we have convergence when /3 < 1; for a = we have stability 
for any /3; and for a < 0, when /3 < 1/2 (convergence for the critical value 
(3=1/2 having been proved above). 



8 



3 Proof of the nonlinear instability result 



As mentioned in the Introduction, a resolvent estimate will be required, we 
prove it first. Then we use an asymptotic expansion of solutions of (NSe) 
around a linearly unstable shear profile Ug, as in [T3]. Like for Theorem [1] 
there will be a slight difference between the positive case (with a > 0) and the 
negative case. An example of linearly unstable shear flow will be given in the 
final paragraph. 



3.1 Preliminary results on the linearised Euler equation 
3.1.1 Linear instability of the Euler equation 

We set a{k) the highest real part of complex numbers A that are unstable eigen- 
values of the linear Euler equation; for each wave number fc, it is finite, and a is 
an even analytic function of fc. In turn, the function a has a maximum ctq > 
(see the study of the Rayleigh equation in [13] , paragraph 4) . We assume that 
this maximum is nondegenerate. 

Fix a wave number fc, and consider the space Vk of divergence-free Fourier 
modes written as u{t,x,y) = v{t,y)e'^'^^ , with v{t,-) G i/" for every s > 0. We 
write divergence-free fields tangent to the boundary by using stream functions: 

uit, X, y) = e"=-(9,vE'(t, y), ~zfcvl/(i, y)) ^ V^(e"=-vl/(i, y)) 

with ^"(^,0) = 0. Note that we have rot u = -A(e*'=^^'(i, y)). For u e Vfc, we 
use the norm ||m||j = \J -I- ||rot . As rot m = (fc^ — dyy)'^, standard 
elliptic regularity (see [H] for example) gives us 

l|V^i(i)IL.(o) = ||V'(e^'^*)||^.(^) ^ llrot u(OIL.(f,) (13) 

Using this, and standard Sobolev inequalities, we get that if w e Vk and u' £ Vk' , 
then u ■ Vm' G Vk+k' and 

\\u-\/u'l<C\\u\\,^Ju'l^^ (14) 
Theorem 3. Let A' > ctq; and ■w{t,x,y) G Vk such that 

\\wmi < ^» (TpTp (15) 

for every I > 2 and for some a > 0. Consider the linearised Euler equation with 
source term 



(ELS) 



dtu + Us ■ Vu -I- u ■ Vug + Vp = w 

div u =0 

u|t=o =uo eVk 

U2\v=0 =Ui eVk 



with |ui(t)| < Cu,(l +t)^'^e^ *, and likewise for dtUi. The initial condition uq 
must be compatible with the boundary condition ui . Then the solution of this 
system satisfies the estimate 

||«(t)IL_2 <C(l+t)""e^'* 
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for t > 0, with u{t) G Vk, and C , which depends on w, Uq, ui, I and k, and is 
locally bounded in the parameter k. 

We show the result for functions with a single Fourier mode, but it extends 
to wave packets 



U{t,x,y) — / ip{k)u{k;t,x,y) dk 
Jr 

with u{k) (z Vk and if a smooth compact-supported function. 

Proof: we adapt the arguments used to prove the resolvent estimate for the 
KdV equation in [3T]. Set 

u{t, X, y) = e"=- ( \ and w{t, x, y) = e^''^w{t, y) 



'ik^{t,y) 

and note that rot u = (fc^ — ^yy)^ ^fc^: where, for fc 7^ 0, the differential 
operator Bu : 7J'+2(R+) n Hl{M+) H\m+) is invertible. 

We first examine ||rot u\\^i. Taking the vorticity of (ELS), we obtain the 
equation 

dtBk-^ + ikiusBk-^ + <'*) = rot w (16) 

We then set 

^^Q{t, y) = *(i, y) - e-^**(0, y) - e~'^^*(t, 0) + e-'^(*+2')*(0, 0) 
with arbitrary ^ > 0, thus ^'o(0, y) = and ^'o(t, 0) = 0. *o solves 
dtBk^o + ik{usBk^o + w"*o) = Fo 

We will not give the detailed expression of Fq, but we point out that integrating 
the hermitian dot-product of Fq by "ifo by parts leads to 

|(Fo|*o)| <c(||/o||i. + ||V,*o||i. 

where fo contains data terms w, uq and ui, and ||Vfc^'o||2,2 is the norm of the 
velocity: 

Vfc*o = (^^-0,9^*0) 

Note that /o = *&i + $26"^*, where $1 contains rot w and terms with y = 
(namely ui and dtui), while $2 contains terms with t = 0. 

Using the Laplace transform in time, Ctg{z, y) — J^°° e~^^g{s, y) ds, and 
writing ^E* = £(4'o and F — CfFo, we turn this differential equation into the 
eigenvalue problem: 

zBk^ + ik{u,Bk^ + <*) = F (17) 

We choose 70 sJctq, A'[, and set z — ^0 + it. 70 being fixed, we abbreviate 
^(z) = ^(t)- When r evolves in M, we get the following estimates: 

Lemma 4. Let I > 0, and ^> solve |_?7| ) with ^'(t, 0) = 0, and F verifying 

mn<c(\\vk^\l^ + \\f\\l?j (18) 
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where f = Ctfo depends on the data. There exists C depending on k, I and 
locally bounded in the parameter k, such that 



Bk^ij) <C||/(r)||^,+. and 



if' 



HI 



<c\\f\\%^. 



Note that the lemma provides estimates on both the vorticity and the ve- 
locity. We prove this lemma in the next sub-paragraph. 
By Parseval's equality, we have 



r+oo r+oo 

/ e-2^°*||Bfc*o||^< dt = / Bk^ir) dr 

Jo J -oo 



W 



— oo 

/+00 
-OO 

< C / e-^^«'\\^^{t)\\%^, dt 

Jo 

r+oo 

+ / e-2(T°+'')*||$2||^,+, dt 
Jo 

r+oo 

< C e-'''^\\\<^,{t)\\%^^ + \\^2\\U.)dt 

Jo 

Replacing $i by $i1[o.t]xr+ for some T > does not affect the solution on 
[0,r] X M+, so iH]) gives us 

pT 2(A'-7o)t 

e-27ot IIBfc^of dt<C dt 



(1 + ^)" 



where C depends on ||u(0)||;, 70, I and fc, and is a locally bounded function of 
k. Noticing that 

(K(OI + |*(o,o)|) 

we get 

T ^ rT 2(A'-7o)t f.2(X'~jo)T 

e-^''°'\\Bkm„i dt<C —dt<C- -— (19) 

II iiH - il + t)°' - (1 + T)« ^ ^ 

(the last inequality is obtained by integrating by parts). 

Note that, using the same procedure, estimate also holds for 



e-'-^^Wmil. + Wdynw) dt < ' ' °. (20) 



10 (1 + T)" 

A quick energy estimate on (fT6|) gives us 

\\B,^m% < c{\\Bk^ml' + Mt)\\% + iirot wm%) 

Multiply this by e^^'^"*, and using and the hypothesis on w, we get 

/ 2(A'-7o)t' 
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Finally, we integrate between and T, and use (|19p and (PU)) to we get the 
desired estimate: ||rot < C(l + t)-"e^'*. 

It only remains to estimate ||w||j;^2. We simply multiply the linearised Euler 
equation (ELS) by u, and get 



ldt\\u\\l..<Ci\\u"' 



L2 -r- l|Vu||i2 + ||w||i2) 



where C depends only on u^. Multiplying by e ^'''o*, integrating in time and 
using ((ni), (Ull) and we get \\u{t)\\^2 < C(l + t)-°e^'*. □ 



3.1.2 Proof of the resolvent estimate (Lemma [4|) 

First note that the estimate is obvious when fc = 0. We must then prove the 
estimate for fc 7^ and then show that the constant C is a locally bounded 
function of fc as fc — > 0. We divide the proof in two parts: with z = 70 + ir and 
T G M, we show the estimate for |r| > M, with M to be chosen independent 
of I, and then show the estimate for |r| < M, using the notion of exponential 
dichotomy (see [5]). 

Step 1: we start with an estimate, and we separate the case |fc| large 
from the case |fc| close to zero. Take |fc| > fco to be chosen later, and the 
imaginary part of the dot-product of equation P7)l and Q = Bk^. The 
elliptic regularity estimate ||^'||^2 < max(l,fc(^^) ||i?/j\E'||^2 then gives us 



C\k\ 



lei 



L2 



with C depending on Ug and fco- Setting ^(t) = |r| — |fc| ||us|lioo — C|fc| 
M so that g{M) = 1, we have, for |t| > M, ||e||^2 < |iF||^2 < 



i and 



Now take |fc| < feg. The imaginary part of the dot-product of ([T7)) and Q 
leads to 



1611^2 <(a|A:| + l)||e||^. 



-\\F"^ 



L2 



L2 



(21) 



using the Young inequality, and that 





2 




fc^ 


< 






L2 



Using the inequality 



(fT5)) . we consider the real part of the dot-product of ([TC)) and 'J: as (u"4'|\E') is 
real, we get 



70 



< 



L2 



C\k\ 



To 
2 



L2 



+ cii/iii2 



Thus, for |fc| < fco such that Cfco 



To 
2 



we have 



L2 



< CI 



with C depending only on 70 and Ug- This allows us to conclude the estimate 
(ED): 

Mr) 11611^2 = (|r|-C.|fc|-l) 11911^2 <C 11/11^2 

and again, we choose |r| > M so that h{M) = 1. Note that M depends affinely 
on |fc|, so we have a uniform bound on ||0(r)||^2 for |fc| < p and |t| > M{p). 
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iJ' estimates are easily obtained by induction: examine the dot-product 
of e?' pT|) by 9' 8, and use Young's inequahty and the induction hypothesis 
lie||^,-i < Cll^^ll^,-! to get f ||5'e||5^, < C||F||^,, with C depending on 
7o, \k\ and I. 

Step 2: we seek a solution to ([T71) as = ^'i + 5*2, where 

-zBk^i = F (22) 
-(z + ^fc^i,)Sfc^'2 + ^fc<*2 = -ifcKBfc*! fcG (23) 

Immediately, we get < :^ and multiplying ([^ by ^i, we get 

^i\\Vk^i\\l. < -{\\'^k^i\\lAt=o + \\F\\l.) 
2 70 

so a ||vl'2|li/i < C ||fcG||^i estimate suffices to prove the result, as 

llfcGII^, < c\k\ + c' ||Vfc*i||^, 

First rewrite (|23p as an ordinary differential system: 

dyU(v)^A{k,z-v)U{v) + H{y) (24) 
with U = (*2, 9,y*2), H ^ (0, kG) and 

/ 1 \ 

A{k,z;y)^\ ^ tku'Jjy) „ ] = A^oik) + B{k, z;y) 

\ ^ z + iku,{y) ^ J 

where Aoo{k) = liuiy^+oa A{k, z;y), and B{k,z;y) = 0{e^^y) for some 77 > 
since u" decays exponentially. The eigenvalues of Aoo{k) are real, so by the 
roughness of exponential dichotomy [5], the system dyU = A{k,z)U has an 
exponential dichotomy on M+; this means that if T{k, z; y, y') is the fundamental 
solution of this last equation with T{k, z;y,y) = I2, there exist a projection 
P{k, z\ y), verifying r(fc, z; y, y')P{k, z; y') = P{k, z; y)T{k, z; y, y'), and positive 
constants C and a, with all of these depending smoothly on {k,z), such that 
for any ^ G C^, 

\Tik,z;y,y')P(k,z;y')^\ < C(fc, z)e-"(^-^') |C| if y > > 
\T(k,z;y,y'){I-P{k,z;y'm < C(fc, z)e«(^-^')|e| if > y > 

For any < p < p' , C{k, z) is uniformly bounded in the set 

^P,P' = {(fc,7o + *T) I |fc| e [p,p'] and |t| < Af} 

but we would like C{k, z) to be bounded in /Co.p' for a certain /?'. 

To get that C(fc, z) is bounded near fc = 0, we follow the proof of the 
persistence of ordinary dichotomies as done in J. Aoo{k) cannot be uniformly 
diagonalised near fc = 0, because the basis of diagonalisation for fc 7^ is 

Instead, we change to the basis 



-Ifcl Ifcl 



ivi{k),V2ik)) 



1 

-|fc| 1 
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vi is an eigenvector of Aoo{k), with eigenvalue — |fc|. Aoo(fc) is therefore trian- 
gular in this basis, and vi spans a space of exponentially decreasing solutions 
of the equation dU = Aoo{k)U. There exists a projection 11 on this space such 
that 



\T^{k;y,y'){I-n{k;y'M < e'l'^Kf-^') |^| si y > y' > 
\T^{k;y,y')U{k;y')^\ < el'^K^"^') |^| si y' > y > 

Consider E{k), the Banach space of functions V such that 



(25) 



E{k) - 

and S the linear mapping defined by 



< +00 



SViy)^ / T^ik;y,y')il -Uik;y'))Bik,z;y')Viy') dy' 



+00 



r^(fc; y, y')n{k; y')B{k, z; y')V{y') dy' 



Note that \B{k,z;y)\ < 6|fc|e"™. Using (gl]), we get that if F e E{k), then 
SV e E{k) with the estimate 



\SV{y)e^''^y\ < \\V\\ 



E(k) 



+00 



b\k\e 



nv 



dy' ) < 2 ll^lU(fe) 



(26) 



when |fc| < p' small enough. So S" is a contracting endomorphism of E(k). 

Let Uoo{k) G E{k) be a solution of dU = Aoo(k)U. By the Duhamel for- 
mula, a bounded solution of the equation dU — AU is a fixed point of the afhne 
transform S — S + Uoo] thanks to Picard's fixed point theorem allows us 
to conclude that such a solution U exists in E{k). Finally, we must uniformly 
bound llC^ll^^j,-). Choosing Uao{k) (e^l''!^, — |fc|e~l'^l^), the decreasing eigen- 
function of Uooi (|26l) gives us the wanted bound. 



The end of the proof is the same as in [JT]; we provide it for completeness. 
Note that for fc ^ 0, / — P{k,z) is the projection on the subspace of solutions 
that go to when y — > -|-oo; let us define Q{k, z) the projection on the subspace 
of solutions such that the first component vanishes at y — 0. As the linearised 
Euler equation does not have an eigenvalue with real part 70, necessarily 

n{i - p{k, z- 0)) n 7^(g(fc, z- o)) = {o} (27) 

and we define a basis of solutions (ei(fc, z), e2(fc, z)), with ei G TZ{I — -P(fc, z; 0)) 
and 62 G TZ(Q{k, z; 0)). Define a new projection 

P'(fc, z) = (ei(fc, z), e2(fc, z)) J (ei(fc, z), e2(fc, z))"i 

and F'(fc, z; y) = T{k, z; y, 0)F'(fc, z), so that 7^(P'(fc, z; j/)) = 7^(Q(fc, z; y)) and 
7?.(/ — P'(fc, z; y)) = 7?.(/ — P{k, z; y)). We also have the estimates 

|T(fc,z;y,y')^'(fc,^;y')ei < C'(fc, z)e-"(f-S'') j^l si y > y' > 
|T(fc,z;y,y')(/-^'(fc,^;y'))ei < C'(fc, z)e"(«-«') 1^1 si y' > y > 

(28) 
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Again, we must check that C"(fc, z) is bounded in A^o.p- To do so, we point out 
that the projections Q{k, z) and / — P(/c, z) can be continued up to A; = 0: we 
have that Q(0, z) is the subspace of solutions to the equation z'^" — with 
^'(0) = 0, and / — P(0,z) is the subspace of bounded solutions to z^" — 
(constants). As the only bounded solution of z^j/" = Q with ^1/(0) = is * = 0, 
(P7)) is true up to A: = 0, so C"(fc, z) is bounded for k near 0. 

Finally, by Duhamel's formula, a bounded solution of is, for fixed [k, z), 

ry r+oo 
U{y) = / ny, y')P'{y')H{y') dy' - / T(y, y'){I - P{y'))H{y') dy' 

Jo Jy 

since the only bounded solution to dyll = A{k, z)U is zero. So, by using ((28)) 
and standard convolution estimates, we have ||J7||^2 < C{k,z) \\H\\^2- To get 
estimates on the norms of the derivatives, just differentiate the equation I 
times, notice that it is the same type of system and conclude by induction. □ 



3.2 Proof of the instability result (Theorem [2]) 

Until the last sub-paragraph, we consider the positive case, which is that of 
the (NSe) equation with the (NCae) Navier boundary condition with a > 0. 
First and foremost, as the initial profile depends on e~'^^^y, we rescale the 
variables {t,x,y) — s^^^{t, x,y), and these new variables will be indifferently 
noted {t,x,y). From now on, the faster variable will be 1" = e^^^'^y. Indeed, 
the system we get after rescaling is 



(NS') : 



div =0 

w^(0,a;,2/) ^ul{x,y) 

(NP): ^^1(^,2^,0) =0 

(NC + 1): \dyu\{t,x,{)) =wf(t,x,0) 



with the initial condition Uq close to a linearly unstable shear flow Us- The 
Navier boundary condition no longer depends on the viscosity, and we have the 
following asymptotic expansion (from 15 ): 

u'{t, x, y) - u7(t, y) + u\t, x, y) + ei/V(t, x, Y) + 0{e'^/^) (29) 

where solves the Euler equation, v!' is the boundary layer and ©(e^/^) is the 
amplitude of the error in L°°([0, T], ^^(n)) for any T > 0. 

Starting with a linearly unstable solution to the Euler equation, we will 
construct an approximate solution to the Navier-Stokes system (NS') above 
that will allow us to prove the instability inequality ([2]). 



3.2.1 Building an approximate solution 

We write =Tl2+ X^^i ^^^U'^ and p"'^ = where TV is large and 

to be chosen later, as an approximate solution of the (NS') system, where each 
will be written as a wave packet. Letting fco be such that cr(fco) = (Tq, and 
setting ip^ as compact-supported smooth functions of k with /cq C supp((^), we 
write 

W{t,x,y)^ [ <f^{k)V^{k;t,y)e''''' dk ^ [ {k)V^ {k;t,x,y) dk 
Jr Js. 
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with V^{k), called terms of the main expansion, being in the previously-defined 
Vk space. We shall write these wave packets more precisely in the next para- 
graph. The equation that the are supposed to solve are Navier-Stokes equa- 
tions linearised around uj: 

dtV^ + III • + ■ VuJ - VsAV^ + VP^ + V^' ■ VV,2 = (30) 

jl +32 =j 

We will not necessarily have div =0, but in total we must have div = 0. 
As we can only solve pop approximately, we write a sub-expansion 

8n-l 

y] = e"/8t,*.8J"+™(t,a;,y)-H£('"+2'/^i;'''^^"+"(i,x,y) (31) 

m=0 

following the expansion (1^^ . and likewise for the pressure - there is a gap of 
e^/'^ between corresponding internal and boundary layers. The first term, y^, 
solves the linearised Navier-Stokes equation above with an error: 

For i > 1, we ask to take the previous error into account, but the fact that 
we once again solve the equation for approximately creates a new error. So, 
solves 

31+32=] 

Fixing a wave number k and writing q = 8jn + m, we must understand the 
equations that w*/'''? = w*/''''(fc) solve to get estimates on these functions. 

We plug the ansatz (PTI) into each equation of the linearised Navier-Stokes 
system, and consider the equations obtained when asking terms of a same order 
e"?/^ to cancel out. We get that the internal layers w*'', which solve only the 
Euler part of the (NS-^e) equation, are solutions of linear Euler systems 

((Eg): dtv'-i + Us -Vv'-i +v''i ■Vus+r'''' + Vp''i =0 
(EL(g)) : \ div v''i = 

[ (Pq): w*'«(t,a;,0)+w^'«-2(t,a;,0) =0 

where we linearise around Ug instead of uj, because we have a better under- 
standing of the linear equation around time-independent profiles, and 

Sn+qi+q2=q 

We also get that the boundary layers each solve a Stokes problem with Navier 
boundary conditions 

{(Sg) : dtv''''^ + v^'i ■ VuS ~ w''''' + r^^i + Vp'''« = 
[Ng) : + ^dyvY" - ~ v'^^-^) (t, x,0) =0 
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the second boundary condition being that v^'"^ should vanish as y — > +00. We 
do not give any detail on r''''^ other than that depends on £~^/*m7 and terms 
u*'"* and w'''"* with m < q. 

In both systems, initial conditions arc chosen to be compatible with the 
boundary conditions, and rapidly decreasing functions of y (like z;(0, 0)e~^); 
therefore V^\t=o is 0{e^). For low values of q, any term with a non-positive 
index is of course ignored, and we must add to the approximate solution a cor- 
rective term £(''^+^)"+^/®w(t, x, Y), to deal with the remainders of the divergence 
equation and the (NP) boundary condition: 

J- a 6,M-1 1/8 Q b.M i (4. n\ f 6.M-1 1/8 b,M\,, n^ 

aw uj = —OxV{ —e' OxV{ and uj2[t,x,0) = {—V2 —e' v^ )(t,a;,0) 

where M = 8nN. Thus the total, u'^p, solves the Navier-Stokes equation with an 
error term and an initial condition e"-close to Ug in H^, but is divergence- free 
and satisfies all the boundary conditions of the (NS') system, although each 
individual does not. 

Now we know the equations that each layer is supposed to solve, we can 
construct the approximate solution by induction. As announced, we start by 
choosing f*''' as an imstablc solution of the linear Eulcr equation (£'1) with 
homogeneous non-penetration condition. In fact, we choose the most unstable 
mode for each k, in other words 

X, y) = e^W*t;(fc; y)e*'=^ 

with SR(A(A;)) = cr(fc). Then, (Di) and the limit condition allow us to choose 
v^'^ = 0, then p^'^ = and iij'" solves a heat equation. 

Having built i)*'"* and v^''^ for m < q, v^''' is the solution of the linear Euler 
system (E(g)). Now we build v^'"^ by taking the equations of the Stokes system 
(S(g)) one after another. First, [Dq] is the equation relative to the divergence- 
free condition: integrating it between Y and -|-oo, we get an expression of the 
normal component of w'''' that goes to as F ^ -l-oo; the rapid-decrease prop- 
erty will follow from the induction below, and nothing guarantees that we have 
V2''^{t,x,0) = 0, so this is where the boundary condition (P,) in the Euler sys- 
tem comes from. Then, notice that thanks to the structure of tZJ, the normal 
part of the Stokes equation {Sq) depends only on v^'' and p'''' (and of course 
previous terms), so we get p^''^ by the same method. Finally, the first component 
of {Sq) is now a heat equation, with a boundary condition given by the Navier 
condition {Nq). 

3.2.2 Estimates on the approximate solution 

The important estimates are those on the first term , or more precisely, those 
on the wave packet of unstable modes of the Euler equation 

u''°(<,x,y)= ( (p^{ky'°{k;t,x,y) dk 
Jm 

with (fi^ supported in / =] — fco — + '7M^o — J?, + where r] is small 

enough so that I C {k \ a{k) > 0}, and such that fco is the only critical point 
of a in ]ko — r],ko + ri[. By Parseval's equality, we have 

lh''°W|li== / vHk)'\\v'''ik;t)\\l, dk^ f ip\kf\\v{k;y)\\l.e^''^''^*dk 
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so, as ib/co is the only nondegenerate maximum of a in each sub-interval of /, 
using the Laplace method, we get ^ f-i/2^2aot ^]-^gj^ ^ _^ +00, so 

we have, given that is bounded near t — 0, 



C 



(32) 



This allows us to estimate the wave packets associated with the other terms 
of U^, which are the u*''? and u''''^ with q < 8n; this is done by induction. Fix 
j < 8n, and suppose that the internal layers u*'' with < q < j satisfy 



, - ,, exp [cto fl + ^) t] 



(33) 



As e'"°"*(l + t)-" +00, we have 



'it)\\H^<C: 



exp [ap (1 + ^) t] 

(1 +i)l/4(l+j78n) 



(34) 



Note that we are studying the stability oiu^{^/et), and that by energy estimates 
on the heat equation with the Robin boundary condition induced by the Navier 
condition, 



rl/8 



< 



1 



rl/8 



so £ 



^^^{u^{y/et) — Us) is bounded in iJ' uniformly for t < e~^^^^. Therefore, 
as the remainder r*^^ in the Euler equation (Ej) verifies as does dtr^'^ , u*'^ 
also satisfies by Theorem [3] 

Now we consider the boundary layer u^'^ . Suppose that u^'', q < j, are 



rapidly decreasing functions in the K-variable and satisfy (I33|) . By construc- 
tion, Wj'"' and p^'^ are rapidly decreasing functions in the y-variable. Then, as 
£~^/^|?I7(t, Y) \ is also uniformly bounded in e for t < £~^/'^^, Y € R+ and £ < £0 
small enough (however there is no guarantee that £~^/*Vu7 verifies such an es- 
timate, so we have to leave it as a term in the equation (Sj)), we get, by using 
the Green function of the heat equation with Neumann boundary conditions 
(see [33]), we get that u'''^ is a rapidly decreasing function of Y satisfying 



< cVTTt 



gO-o(l+0-l)/8n)t 
(1 +i)l/4[l + 0'-l)/8n] 



exp [ap (1 + ^) t] 

(1 +i)l/4(l-hi/8n) 



Overall, we have 



k"f^^(i)| 



L2 



£ 



n+j/8 



1^ fl _Lt)l/4(l+j78«) 



exp 



<yo 1 



8n 



(35) 



But we will be interested in times t = Tp — t, with Tp such that 
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The important manoeuvre here is to write i = — r, so that 



1 + 



y 1 + - T 

with Kf^{T) g]1,2] for e small enough, and with < r < tq where tq does not 
depend on e. As Tg ~ cln(e~^), we can choose e < eg so that Tq — tq > 0. tq, 
and therefore will be chosen at the end of the proof in the next paragraph. 

We show how this manoeuvre works in detail. Writing i = — r, ([55)) 
becomes 

8n-l 

||e"f/i(TJ - t)||^, < (1 + - TY'^C^Kl{T)e-^^'^ ^ (/^o^(r)e--«^)-'"/8" 

< (1 + - r)i/H6nCii^g(T)e-'^°^ 
Thus, by returning to the t- variable, for t <Tq, 

||[/iW||^,<(l+0^/^^^ (36) 

Note that this estimate also means that ||f^"'^(i)||i2 behaves like ||u''°(t)||^2- 

In Theorem [2l we want an Z/°° under-estimate, so we need to work more 
in order to get an under-estimate of a local l? norm. If a = 3(A'(A;o)) and 
^ = -cr"(fco) > 0, then 



|M^-"(t,x,y)|^e-''* 



ip^{k)v{k,y) exp{i{x + at){k - ko) - (3t{k - fco)^) dk 



But, writing u = k — ko, we have (by factorising) 

exp(i(a; + at)u — ptu ) du ~ — -= exp 



with c complex, so there exists C such that 



^fin M ^ ^e'^"* „ (-{x + atf 



W'^{t,x,y)\ > ^— — exp 



Integrating on the bounded domain 



nA{t) = {(x, y) \ y< Aa,nd\x + at\ < A\/l + t} 

for A large enough, and using the fact that u''*^ is the dominant term in U^, we 
have 

II^^WlL^(n.,))>^Ul+t)V4-^ (37) 



As the measure of ftA{t) is A^^/1 + 1, this gives us an under-estimate of the 
norm ||C/H^)|Loo- In fact, ||C/^(i)|Loo ^ C{1 + t)~^^^e''°K 

Now we can get estimates on for j > 1. Remember that behaves like 
v^'°, so induction using Theorem [3] gives us 

\\vHk;t)\\^^<c,e^^ot (38) 
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Indeed, for j > 1, the internal part of solves a linearised Euler equation with 
a source whose predominant terms are V^^ ■ ^V^^ with ji + j2 — j, which verify 
([55)1 by (HH) , and as above, the boundary layers are rapidly decreasing functions 
of Y satisfying p3[) . with q > 8n. Now we specify the structure of U^: we write 



V-' {ki, ■ ■ ■ ,kj]t,x,y) dki ■ ■ ■ dkj 



with \V^{ki, ■ ■ ■ ,kj;t,x,y)\ < C exp[(cr(fci) + • • • + (T{kj))t], which gives us 



exp[(f7(fci) + • • • + a{kj))t] dki ■ ■ ■ dkj 



\U'{t)L^^ < C 



< C l^J exp{aot ~ I3{k ~ koYt) dk ] < 

thanks to the second-order Taylor inequality 

a(fc) < do - - fco)' 
For an estimate, we rewrite as 



(1 + ty/^ 



(39) 



(40) 



Wit,x,y) 



T/^(fei,-..,fc,;i,2/)e^^'=- 



Ijkejl Jki-\ \-kj=jk 

with jl = I + ■ ■ ■ + I. By Parseval's equality, we have 



\U\t)\ 



jkeji 



ki-\ \-kj—jk 



V\k,,---,k,;t) 



L2 



]k£jl 



Using (|40|) . and noticing that 



m— 1 



m— 1 



we have 



N{tf < Ce(^''^''-^''5('=-'=«)')* / exp(-/3 V(fc™ - k)h) 



52 krrt=jk 

As fcj is imposed by the other j ~ 1 fcm-variables, we integrate j — 1 gaussian 
functions, so 

(Jg2jaot-2jl3{k-ko)^t 



\uHt)\\ 



L2 



< 



p 



i-1 



dk 



Integrating this final gaussian function, and taking into account the boundedness 
of II ^^"'(0 11^2 near t = 0, we have 



(1 + ty/^ 



(41) 



This easily gives us a final estimate on the remainder of the equation that 
u°-P solves: 

dtu^P + u^P ■ Vu^P - ^/lAu^P + Vp'^P = R 
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The remainder R contains terms of the type ■ VU^ with j + 1 > N, so by 

using (|14p . we have the estimate 



2N 



Replacing t by Tg — t, we get 



j=N+l ^ ' 



2N 

and so, for t e [T^ - tq, T^], 

< + 0^/^e(^+i)-' (42) 

Cfl can be adjusted so that the inequahty is vahd for t £ [0,Tq]. 
3.2.3 Proof of the instability 

For this paragraph, we denote S = e", and we choose N such that for t > 0, 
2Nao > 2 ||Vu7(t)||^^ + 1. Let be an exact Leray solution to the (NS') 
system with initial condition ^"^(0). We have ||m^(0) — Us|li2 < e"- The aim is 
to get a lower bound of \\u'^ — uJH^^o, so we will look for a local lower bound. 
Let us first estimate — ti7|lL2(j2^(t))- using (1411) and ((37)) . we have 

Writing i = — r, we get 

with Cg = -/V2^ maxi<j<Ar Cj. There exists ri > such that 

for r > Ti. So, choosing tq > ti, we have, for Tq — tq < t < Tq — ti, 

\\u^nt)~^mmnAt)) > (^+^y^ \^i%i/2 ^"' (43) 

Note that ti does not depend on 6. 

Now we estimate ||w(t)||^2(-s-j'), where w = u'^ — u"-^ and q = — v°'^ ■ The 
pair (w, g) solves the equation 

dtw + w • Vw - VeAw + Vg = --u^^ • Vw - w • Vii"'' + R 
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We compute the dot-product of this equation and w; integrating by parts, 
using the divergence-free condition on the nonhnear term, the (NC -I- 1) and 
(NP) boundary conditions on the Laplacian term and Young's inequahty on 
the remainder, we get the energy estimate 



dt \\w\\l. -f 2^e\\Vw\\l.+A^e f \wi\' < (2 || V^.'^^'H^^ + p) \\w\\l. + -\\R\ 

J an P 



2 

L2 



with < p < 1/3. The L°° estimate on Vm"^ is crucial here: as each can 
be written — (t, x, y) + e^^^U^'^ (i, x, F), the first order derivatives of 
verify dm-type estimates, since e^'^V iU^'^ {t,x,£-^^^y)) = (VC/*^^)(i, cc, F), so 



\Vu^^{t)\\^^<\\u;{t)\\^^+Y.7i 



Using the t = — t manoeuvre, the sum on the right is smaller than p for 
t < Tg — T2, with T2 > Ti independent from J, and with (l42t and the choice of 
N we made at the beginning of the paragragh, the energy estimate becomes 

Therefore, 

by applying the following lemma, which is obtained using a variant of Gronwall's 
lemma from |26j and integration by parts: 

Lemma 5. Let ip be a function such that 

dMt) < Xip{t) + C- 



(1 + i)" 

for every t > 0, and for parameters fi > X > and a > 1, then 

ip{t) < C- 

^ - (1 + 
with C depending on X, p and C . 

With (pS)) and (HU, we can now conclude: 

\\u'{t) ~ > Wu^^t) - Il7(i)IL2(a^(,)) - lkWllL2(a) 



> {l + t) 



1/4 



_2(l+<)l/2 (l+;)(A'+l)/2 

Writing t = Tq — t for the last time, we have 



2(l + i)i/2 (i + i)(N+i)/2 -4(i + t)i/2 
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for t < Tg — T3, with T3 > T2 independent from S. Setting T"^ = Tq ~ T, with 
T3 < T < To fixed (here we can finally fix tq and Eq), we have 

\\u^T-) - II7(T^)IL.(o,(T.)) > (1 + T-)i/4^e--''^ (1 + T-^/^So 

where Sq does not depend on e. Thus, ||M^(r^) ~'U^{T'^)\\^^ > C6o- 

Also, as ^ L°°, we have ||m^(T'^) - u7(T^)||^, > CSq. Remember that 
we are still using the fast variables, so, returning to the original scale of time 
and space (again without changing notation), we get ([2]) and 

In the original scale of time, = C'(ny^ln(e^^)). □ 



3.2.4 The negative case 

The proof in the negative case is exactly the same up until the energy estimate 
of w. Here, the Laplacian term becomes 



Aw ■ w — llVw 



n 



2(o) -2|a|||u;i||^2(an) 



The second term cannot be ignored, so we use the same idea as in paragraph 2: 
the trace theorem, Sobolev interpolation and the Young inequality allow us to 
get 

4|a|Ve|lw|1^2(af^) < 2^/e ||Vw||^2(o) + "^CaV^ \\w\\l2(n) 

where Ca does not depend on e. Choosing N so that Nao > || Vu7||^oo +Ca + ^, 
everything works as above. □ 



3.3 An example of unstable profile 

Examples of piecewise-linear flows that are linearly unstable for the Euler equa- 
tion are given in [7 and 13 , the latter stating that close-enough regularisations 
of these are smooth unstable shear flows. We wish to provide an explicit exam- 
ple of smooth linearly unstable profile, which is the object of Proposition [6] 
below. From this example, we easily deduce the expression of a smooth unstable 
profile that fits Theorem [2j 

Proposition 6. For every S > and C € R, the shear flow {us, 0) with 

us{y) = tanh(2/ - S) + ( 
is a linearly unstable for the Euler equation. 

Proof: writing A — —ike for 3?(A) > and k ^ 0, and taking the curl 
of equation (EL), the Euler equation linearised around Ug, we obtain a linear 
second-order differential equation for ^P: the Rayleigh equation 

(R(c, k)) : {us - c)(dly - - <vl/ ^ 
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with the conditions ^"(0) — and hniy_y+oo ^{y) = 0. Our problem is now 
finding c in the complex upper-half-plane, and real numbers k such that (R(c, k)) 
has a solution in H^{R+) n ij2(R+). 

It is well known that if the Rayleigh equation has an unstable solution, then 
Us must have an inflection point in R"*"* (Rayleigh's theorem, ^7\), and us has 
exactly one inflection point, yo = S, with inflection value uq = us{yo) = C- Also, 
the function 

Us(y)-uo tanh(j/-d) 

has a limit when y —> yo, and is a positive continuous function on R+, vanishing 
at infinity. Thus, choosing c = wq, we can divide (R(c, fc)) by Ug — uq (which is 
not usually possible when c is in the range of Us), and we have a Sturm-Liouville 
problem 

^^'\y)-Ksiy)^{y) = -e^{y) 

in which the square of the wave number k intervenes as an eigenvalue of the 
operator Sg — —d^ — Kg. We shall therefore use the following result by Z.Lin 
(Theorem 1.5 in P]): 

Theorem 7. Let U be a profile which has a limit I as y +oo, and an 
inflection point yo such that, writing uq = U{yo), the function K{y) = ^^-^^^^^^ 
is a positive continuous function which goes to zero as y ^ +oo. U{y) can 
take the value I only a finite number of times. If the lowest eigenvalue of the 
operator S — — K, defined on the Sobolev space HqCR^) n iJ^(R+), is 
strictly negative, then, letting —pi? be the lowest eigenvalue, for every k £ (0, pL), 
there exists c{k) such that the Rayleigh equation (R(c(fc), fc)) has an unstable 
solution. 

The result is first shown on a finite interval; one proves that eigenfunctions 
of S with negative eigenvalues are limits of unstable solutions, and this is done 
by using the Picard fixed point theorem. A compactness argument is then used 
to get the result on the half-line. We detail the proof no further. 

In the case of us, it remains to show that Ss — —d^ — Ks has a negative 
eigenvalue. We hope to use a standard Sturm-Liouville argument, and notice 
that ifsiy) = tanh(y — (5) is a solution to the equation 

-^'{y)~Ks{yMy)^0 (45) 

that has exactly one zero in ]0, -|-oo[. Since Ks decreases exponentially to 0, the 
multiplication by Ks is a compact perturbation of ~d^, so by Weyl's theorem 
(see [IH]), the essential spectrum of Ss is M"*". But that means trouble, because 
the in corresponds to the lowest point in the essential spectrum, and 
the one zero of ips means that Ss has between and 3 negative eigenvalues, 
according to the classical Sturm-Liouville analysis shown in "W , chapter XIII. 

So we must use a different tool to determine that Ss has a strictly negative 
eigenvalue. Fix S and let Qs be the quadratic form associated with Ss: for 
u e H^iM.+ ), we define 

Qsiu)^ / \u'iy)\^-K{y)\uiy)\^ dy 
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Let Vri{y) = tanh(y — 77) for rj > 0, and x be a [0,2]-supported smooth function, 
with x{y) = 1 foi' 2/^1- any n G N*, we define 



wJJ is equal to v^j on [(5, n], continuous on R+, and is in iJg(M+), and since v'^ 
and \/ KsVri are square-integrable on K"*", we have 



By integrating by parts, we have Q{5) — 0, and Q is a differentiable function of 
77, so we look at dj^Q{5). We have 



Thus, since vs{y), v'g{y) and Ks{y) are positive for y > 5, drjQ{S) > 0. So there 
exists 770 such that, for 77 g]7;o, S[, Qi^i]) < 0, and therefore, for a given 77 in that 
interval, there exists n large enough such that verifies Qs{wl^) < 0. 

We have proved that the lowest point of the spectrum of Ss is negative, and 
Us verifies all the hypotheses of Theorem [3 Proposition [6] is proved. □ 

Remark: Ss has in fact exactly one negative eigenvalue. Remember that 
'^s{y) = tanh(7/ — 5) is a bounded continuous solution of —Lp" — Ks^p = 
that vanishes only at 7/0 = 5, and that Qs{u) = J^°° \u'\^ — Ks\u\'^ is the 
quadratic form on iJQ(R+) associated with the operator Ss- We shall show 
that Qs > on Fs = {u G -ffo(R"*") I "(^) = 0}. A function u e Fs can 
be written as u = ipsv, with v £ Hq{M.'^). Let us take v g Cg^(M"'"), with 
supp(7;) C [h, yo — h]U [7/0 + +00 [ for a certain h > 0. Replacing u by tpsv in 
Q{u), we get 

— 7i" — Ksu = —(fisv — 2ip'sv' — ifisv" — Kipsv = —2ip'gv' — ipsv" 

so, using Qs{u) = J^°° Ssu ■ u, we have Qsi'fisv) = Jq°° -2ipsip'svv' - Lp\vv" . 
Integrating by parts, we get the factorisation Qs{fsv) = f^v'^, ie 



as Cq'{{0,6) U {5, +oo[) is dense in Fs- We now have Qs > on Fs, which is a 
hyperplane of Hq{R'^), thus Qs is negative only on a subspace of dimension 1: 
Ss has only one negative eigenvalue. 

We can now give an example of unstable shear flow for Theorem (2) it is one 
of the linearly unstable profiles of the Euler equation, us{y) = tanh{y — <5) + C 
with C and S to be chosen such that it satisfies the rescaled Navier condition. 







for all u e Fs 



(46) 



-^^^^(O) = aus{0) 



(47) 
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This example will cover all cases except a — 0, firstly because u'g never vanishes, 
and secondly because when a = 0, instability does not occur (see [T]). (|T7)) 
is an equation of the variable X ~ tanh((5) with parameter ^; precisely it is 
^X^ — aX + aC — 5 = 0; the discriminant is A = — 2a( + 1. We then have 
two cases: 

• for a ^ Z, we choose ^ so that A = i?[a]^, therefore one of the solutions 
of the polynomial equation is Xq = a — E[a\ g]0, 1[; 

• for a e Z*, we choose ^ such that A — {a — i)^, so Xq — i. 
Setting tanh^^(Xo) Sq, usg is provides the wanted example; explicitly 

^ , ^ f tanh(y- tanh-i (a- i?[a])) +f - Mid!^ if „ ^ z 
' [ tanh(y-tanh"^(i)) +^ + |r if a e Z* 

are flows that fit Theorem [2] 
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